Model order reduction for left ventricular mechanics via congruency training

Computational models of the cardiovascular system and specifically heart function are currently being investigated as analytic tools to assist medical practice and clinical trials. To achieve clinical utility, models should be able to assimilate the diagnostic multi-modality data available for each patient and generate consistent representations of the underlying cardiovascular physiology. While finite element models of the heart can naturally account for patient-specific anatomies reconstructed from medical images, optimizing the many other parameters driving simulated cardiac functions is challenging due to computational complexity. With the goal of streamlining parameter adaptation, in this paper we present a novel, multifidelity strategy for model order reduction of 3-D finite element models of ventricular mechanics. Our approach is centered around well established findings on the similarity between contraction of an isolated muscle and the whole ventricle. Specifically, we demonstrate that simple linear transformations between sarcomere strain (tension) and ventricular volume (pressure) are sufficient to reproduce global pressure-volume outputs of 3-D finite element models even by a reduced model with just a single myocyte unit. We further develop a procedure for congruency training of a surrogate low-order model from multiscale finite elements, and we construct an example of parameter optimization based on medical images. We discuss how the presented approach might be employed to process large datasets of medical images as well as databases of echocardiographic reports, paving the way towards application of heart mechanics models in the clinical practice. © 2020 Di Achille et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.19-14- 00134Russell Sage Foundation, RSFSK and OS were funded by RSF (http:// www.rscf.ru/en/) as described below. Part of this work was carried out within the framework of the IIF UrB RAS government assignment and was partially supported by the UrFU Competitiveness Enhancement Program (agreement 02. A03.21.0006) as well as the RSF grant (No. 19-14- 00134). The Uran supercomputer at IMM UrB RAS was used for part of the model calculations. IBM provided support in the form of salaries for authors PA, JP, JK and VG but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of these authors are articulated in the "author contributions" section

Meshless electrophysiological modeling of cardiac resynchronization therapy—benchmark analysis with finite-element methods in experimental data

Computational models of cardiac electrophysiology are promising tools for reducing the rates of non-response patients suitable for cardiac resynchronization therapy (CRT) by optimizing electrode placement. The majority of computational models in the literature are mesh-based, primarily using the finite element method (FEM). The generation of patient-specific cardiac meshes has traditionally been a tedious task requiring manual intervention and hindering the modeling of a large number of cases. Meshless models can be a valid alternative due to their mesh quality independence. The organization of challenges such as the CRT-EPiggy19, providing unique experimental data as open access, enables benchmarking analysis of different cardiac computational modeling solutions with quantitative metrics. We present a benchmark analysis of a meshless-based method with finite-element methods for the prediction of cardiac electrical patterns in CRT, based on a subset of the CRT-EPiggy19 dataset. A data assimilation strategy was designed to personalize the most relevant parameters of the electrophysiological simulations and identify the optimal CRT lead configuration. The simulation results obtained with the meshless model were equivalent to FEM, with the most relevant aspect for accurate CRT predictions being the parameter personalization strategy (e.g., regional conduction velocity distribution, including the Purkinje system and CRT lead distribution). © 2022 by the authors. Licensee MDPI, Basel, Switzerland

APPLICATION OF NEURAL NETWORKS FOR PROCESSING OF HEART MEDICAL IMAGES

Deep learning has become the most widely used approach for cardiac image segmentation and identifying predictors of its pathologies in recent years. To solve this problem, a modifi-cation of the U-net architecture was used. The segmentation accuracy was about 80%

Machine Learning Prediction of Cardiac Resynchronisation Therapy Response From Combination of Clinical and Model-Driven Data

Background: Up to 30–50% of chronic heart failure patients who underwent cardiac resynchronization therapy (CRT) do not respond to the treatment. Therefore, patient stratification for CRT and optimization of CRT device settings remain a challenge. Objective: The main goal of our study is to develop a predictive model of CRT outcome using a combination of clinical data recorded in patients before CRT and simulations of the response to biventricular (BiV) pacing in personalized computational models of the cardiac electrophysiology. Materials and Methods: Retrospective data from 57 patients who underwent CRT device implantation was utilized. Positive response to CRT was defined by a 10% increase in the left ventricular ejection fraction in a year after implantation. For each patient, an anatomical model of the heart and torso was reconstructed from MRI and CT images and tailored to ECG recorded in the participant. The models were used to compute ventricular activation time, ECG duration and electrical dyssynchrony indices during intrinsic rhythm and BiV pacing from the sites of implanted leads. For building a predictive model of CRT response, we used clinical data recorded before CRT device implantation together with model-derived biomarkers of ventricular excitation in the left bundle branch block mode of activation and under BiV stimulation. Several Machine Learning (ML) classifiers and feature selection algorithms were tested on the hybrid dataset, and the quality of predictors was assessed using the area under receiver operating curve (ROC AUC). The classifiers on the hybrid data were compared with ML models built on clinical data only. Results: The best ML classifier utilizing a hybrid set of clinical and model-driven data demonstrated ROC AUC of 0.82, an accuracy of 0.82, sensitivity of 0.85, and specificity of 0.78, improving quality over that of ML predictors built on clinical data from much larger datasets by more than 0.1. Distance from the LV pacing site to the post-infarction zone and ventricular activation characteristics under BiV pacing were shown as the most relevant model-driven features for CRT response classification. Conclusion: Our results suggest that combination of clinical and model-driven data increases the accuracy of classification models for CRT outcomes. Copyright © 2021 Khamzin, Dokuchaev, Bazhutina, Chumarnaya, Zubarev, Lyubimtseva, Lebedeva, Lebedev, Gurev and Solovyova.This work was supported by Russian Science Foundation grant no. 19-14-00134

Generative Adversarial Networks for Construction of Virtual Populations of Mechanistic Models: Simulations to Study Omecamtiv Mecarbil Action

Biophysical models are increasingly used to gain mechanistic insights by fitting and reproducing experimental and clinical data. The inherent variability in the recorded datasets, however, presents a key challenge. In this study, we present a novel approach, which integrates mechanistic modeling and machine learning to analyze in vitro cardiac mechanics data and solve the inverse problem of model parameter inference. We designed a novel generative adversarial network (GAN) and employed it to construct virtual populations of cardiac ventricular myocyte models in order to study the action of Omecamtiv Mecarbil (OM), a positive cardiac inotrope. Populations of models were calibrated from mechanically unloaded myocyte shortening recordings obtained in experiments on rat myocytes in the presence and absence of OM. The GAN was able to infer model parameters while incorporating prior information about which model parameters OM targets. The generated populations of models reproduced variations in myocyte contraction recorded during in vitro experiments and provided improved understanding of OM’s mechanism of action. Inverse mapping of the experimental data using our approach suggests a novel action of OM, whereby it modifies interactions between myosin and tropomyosin proteins. To validate our approach, the inferred model parameters were used to replicate other in vitro experimental protocols, such as skinned preparations demonstrating an increase in calcium sensitivity and a decrease in the Hill coefficient of the force–calcium (F–Ca) curve under OM action. Our approach thereby facilitated the identification of the mechanistic underpinnings of experimental observations and the exploration of different hypotheses regarding variability in this complex biological system. © 2021, The Author(s).This work was partially supported by the EU’s Horizon 2020 research and innovation program under the Marie Sklodowska-Curie (g.a. 764738) state Program (No. AAAA-A19-119070190064-4) and the research grant from RFBR (No. 19-31-90089)

ELECTROMECHANICAL MODEL OF THE LEFT VENTRICLE OF THE HUMAN HEART

We developed a mathematical model of the left ventricle electromechanical function of the human heart. We examine one full cycle of the left ventricle contraction with various parameters of the speed of the wave of excitation and spin of the myocardial fibers in the ventricle wall

Thermodynamics of the interacting Fermi-system in the Static Fluctuation Approximation

We suggest a new method of calculation of the equilibrium correlation functions of an arbitrary order for the interacting Fermi-gas model in the frame of the static fluctuation approximation (SFA) method. This method based only on the single and controllable approximation allows to obtain the so-called far-distance equations (FDEs). These equations connecting the quantum states of a Fermi particle with variables of the local field operator contains all necessary information related to calculation of the desired correlation functions and basic thermodynamic parameters of the many-body system considered. The basic expressions for the mean energy and heat capacity for electron gas at low temperatures in the limit of high density were obtained. All expressions are given in the units of r_s,where r_s determines the ratio of a mean distance between electrons to the Bohr radius a_0. In these expressions we calculated the terms of the order r_s and r_s^2, correspondingly. It was shown also that the SFA allows to find the terms related with high orders of the decomposition with respect to the parameter r_s.Comment: 22 pages, 5 figure

Oscillations of the inertia moment of a finite fermi system in the cranking model framework

Section III. Theory of Atomic Nucleus and Fundamental Interaction

Superfluidity of heated Fermi systems in static fluctuation approximation

Section III. Theory of Atomic Nucleus and Fundamental Interaction